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Purpose: Monocyte chemoattractant proteins (MCPs) are important cytokines that involved in cellular activation and releasing of inflammatoy me- 
diators by basophils and eosinophils in allergic disease. Some MCPgene variants implicate in asthma and monoclonal antibody for MCP-3 blocks 
allergic inflammations in the patients with asthma. Detection of interactions between gene and environment or between genes for complex disease 
such as asthma is important. We searched for an evidence of genetic effect of single nucleotide polymorphisms (SNPs) of MCPgenes as well as 
gene - gene interactions involved in asthma. Methods: Four hundreds asthmatics and four hundreds normal controls were enrolled. Asthma was 
defined as a positive bronchodilator response or positive methacholine provocation test with compatible clinical symptoms. Seven MCPgene SNPs 
(2 SNPs in MCP-1, 1 in MCP-2, and 4 in MCPS) were included. Association analyses between SNP and asthma, and the tests for gene - gene inter- 
action were performed. Results: Strong linkage disequilibria were found among 7 MCPgene polymorphisms. There was no SNP that showed a 
significant association with asthma among 7 SNPs of 3 MCPgenes. No haplotype was associated with asthma, either. The combination of MCP1- 
2518G>A, MCP2+46A> C, and MCP3+5630 T was the best predictive model for asthma as compared to the control in tests for gene - gene in- 
teraction. The MCP1-2518G>A and MCP2+46A>C was the second best predictive combination and this had the highest synergistic interaction ef- 
fect on the subject's status than any other combination of polymorphisms. Complete linkages were not associated with the gene - gene interactions 
models. Conclusions: MCPgene polymorphisms probably interact with each other; thus, these findings may help in developing a possible genet- 
ic marker to predict asthma. 
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INTRODUCTION 

Asthma is a common respiratory disease characterized by re- 
versible airway obstruction, bronchial hyperresponsiveness, 
and respiratory symptoms resulting from chronic airway in- 
flammation. Recently, it has been focusing on heterogeneity of 
asthma, which is asthma consists of multiple phenotypes. 1 Var- 
ious factors, including cytokines, genes, and other environmen- 
tal exposures, are involved in the development of the disease, 
and allergic inflammation, among many factors, is still impor- 
tant in the pathogenesis of asthma. 2 

Eosinophils are involved in allergic inflammation in asthma 
via several processes, including adhesion to endothelial cells, 
transendothelial migration, and chemotaxis through monocyte 
chemoattractant proteins (MCPs). 3 ' 4 The three MCPs (MCP-1, 
-2, and -3; MIM# 158105, MIM# 602283, and MIM# 158106, re- 
spectively) share 51%-78% amino acid homology and activate 
monocytes via the chemotactic cytokine (CC) P chemokine re- 
ceptor (CCR) 2. MCP-3 activates both CCR1 and CCR3, while 



MCP-2 activates CCR5. 3,4 

MCP-1, -2, and -3 mRNAs are expressed in airway smooth 
muscle cells following stimulation with the cytokines interleu- 
kin-1 beta (IL-1|3), tumor necrosis factor alpha (TNF-tx), and in- 
terferon gamma (IFN-y). Moreover, MCP- 1 mRNA expression is 
associated with an increase in the release of MCP-1 and MCP-2 
proteins. 5 The level of MCP-1 is higher in the bronchoalveolar 
lavage fluid of subject with asthma than in non-asthmatics. 6 
The mean serum MCP-1 level is significantly higher in asthma 
than in control, 7 and this higher level of serum MCP-1 is observed 
during the asymptomatic phase and declines after treatment. 7 
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Furthermore, MCP-1 might play an important role in allergic 
responses via induction of mast cell activation and release of 
leukotriene C4 into the airway, which directly induces airway 
hyperresponsiveness. 8 MCP-3 may play an important role in 
macrophage recruitment in children with upper respiratory vi- 
ral infections. 3 RANTES and MCP-3 mRNA was significantly el- 
evated in bronchial mucosa of atopic asthma, and associated 
with significantly increased numbers of bronchial mucosal eo- 
sinophils. 9 The expressions of MCP-3 in bronchial mucosa were 
also observed together with eosinophilic cytokines such as IL-5 
GM-CSF, and IL-13 in both atopic and nonatopic asthma. 10 

All three human MCP genes are located on chromosome 17q 
11.2-ql2. The frequency of the MCP1-25180A G polymor- 
phism in the gene regulatory region of MCP-1 is significantly 
higher in asthma and correlated with increased eosinophil lev- 
els in Caucasian asthmatic children. 11 On the other hand, there 
is lack of study for the association between MCP-2 and MCP-3 
gene polymorphisms with asthma. 

Because asthma is a complex and heterogeneous respiratory 
disease, the factors involved in the development of asthma, in- 
cluding cytokines, genes, and other environmental exposures, 
likely interact. Predisposition to complex disease imparted by 
individual genes may act independently or interact with other 
genes that exist in the same biological pathway to produce vari- 
able effects. 12 MCP genes may contribute to the development of 
asthma or other allergic phenotypes through receptor-ligand 
interactions. MCP gene structure and function have been char- 
acterized, in addition to sequencing; however, littie is known 
about interactions among MCP genes in asthma. In this study, 
we evaluated the association between the polymorphisms of 
MCP genes and asthma, and further gene - gene interactions 
among MCP gene polymorphisms in patients with asthma. 

MATERIALS AND METHODS 
Subjects 

Four hundred asthmatics and Four hundred normal controls 
were enrolled in this study. The institutional review board of the 
Soonchunhyang University Bucheon Hospitals approved the 
study. All asthmatics had current symptoms, such as wheezing, 
dyspnea, or cough, and met the criteria for asthma as deter- 
mined by the American Thoracic Society. 13 Each patient showed 
airway reversibility as documented by an inhalant bronchodila- 
tor-induced improvement of forced expiratory (FEV1) of > 15% 14 
and/or airway hyperresponsiveness as provocative concentra- 
tion of methacholine required to cause a 20% decrease in FEV1 
(PC20) of <8 mg/mL. 14 Normal subjects were recruited from the 
general population, based on the following criteria: a negative 
screening questionnaire for respiratory symptoms, 15 a predicted 
FEV1 >75%, PC20 >8 mg/mL, total serum IgE <300 IU/mL, 
and a normal chest X-ray. The clinical characteristics of all asth- 
matic and control subjects are presented in Table 1. 



Table 1. Clinical manifestations of asthma and control subjects 





Asthma 

no Li 1 1 1 1 a 


Pnntrnl 

uUI 1 LI Ul 


Pv/ah ip 


No. of subjects 


400 


400 




Age, yr : Median (Range) 


46(8-68) 


29(8-70) 


<0.001 


Sex (Male/Female) 


166/234 


180/220 


0.19 


Smoker (%) 


26.6 


41.6 


0.001 


FVC% (pred) 


88.3 ±1.3 


89.3+1.1 


0.29 


FEV1 % (pred) 


75.7 ±1.2 


96.5±0.7 


<0.001 


% changes of FEV1 by 
bronchodilator 


12.9+0.9 


2.8±0.4 


<0.001 


PC20 (mg/mL) 


2.9±0.2 


24.8 ±0.3 


<0.001 


Total IgE (IU/mL) 


441 .6 ±40.1 


127.6+9.3 


<0.001 



Mean±SEM; Chi-square tests were used for comparisons of categorical vari- 
ables, Student's f-test for continuous variables. 



Genotyping by single base extension and electrophoresis 

Single-base extension reaction was used to genotype poly- 
morphic sites. SNaPshot ddNTP Primer Extension Kit (Applied 
Biosystems, Foster City, CA, USA) was used for primer exten- 
sion reactions according to the manufacturer's recommenda- 
tions. The primer extension reactions were cleaned by the addi- 
tion of 1 U of shrimp alkaline phosphatase and incubation at 
37°C for 1 hour, followed by 15 minutes at 72°C for enzyme in- 
activation. DNA samples, containing extension products, and 
GeneScan 120 Liz size standard solution were added to Hi-Di 
formamide (Applied Biosystems), as recommended by the 
manufacturer. The mixture was incubated at 95°C for 5 minutes, 
cooled on ice for 5 minutes, and electrophoresed using the ABI 
Prism 3100 Genetic Analyzer. Data were analyzed using ABI 
Prism GeneScan 3.5.1 and Genotyper 3.6 software (Applied 
Biosystems). 

Hardy-Weinberg equilibrium, linkage disequilibrium, 
haplotypes, and association analyses between asthma and 
normal control. 

A total of 7 SNPs in MCP-1, -2, and -3 genes (2 SNPs in MCP-1 
gene, 1 in MCP-2, and 4 in MCP-3) were included in the analy- 
sis. Testing for Hardy-Weinberg equilibrium 16 and the calcula- 
tion of D' for identification of linkage disequilibrium were per- 
formed using PDA l.O. 17 Haplotypes were inferred using EM-al- 
gorithm (PHASE 2.0.2) 18 in 2 ways; with separate 3 genes for each 
and all at once. 

Continuous variables, such as serum total IgE levels, and PC20 
methacholine levels were log transformed to approximate nor- 
mal distributions. Chi-square tests were used for comparisons 
of categorical variables, Student's Mest for continuous variables. 
The differences in the log-transformed values among genotypes 
were examined using a logistic regression model. All single SNPs 
and haplotypes (frequency >0.05) were analyzed for their asso- 
ciation with asthma, using logistic regression models with ad- 
justment of age and sex. Statistical significance was defined at 
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the standard 5% level. These analyses were performed using 
software R (www.r-project.org). 

Tests for gene - gene interactions 

To find gene - gene interaction model, we used the muMfac- 
tor dimensionality reduction (MDR) analysis (MDR, v. 10.0). 
Cross-validation is a measure of the number of times a particu- 
lar set of loci is identified in each of the subjects. The threshold 
ratio is defined as the ratio of the number of affected individu- 
als to that of unaffected individuals. The subjects are at high risk 
only when the threshold ratio exceeds 1.0. When the final best 
model was selected, a model for high- and low-risk genotype 
combinations was formed using an adjusted threshold that was 
equal to the ratio of cases and controls in a model that maximiz- 
es the cross-validation consistency (CVC). 19,20 Accuracy is defined 
as the proportion of subjects that are grouped correcuy accord- 
ing to their case-control status. When the CVC was maximal for 
1 model and accuracy was maximal for another, the model with 
fewest loci/factors was selected. 19 21 To test the universal null 
hypothesis of no association of the independent and additive 
main effects, we further used permutation tests with software 
MDRpt. 22 To estimate the contribution of associated genotype 
combination, we calculated the odds ratio of each genotype 
combination in best model using odds ratio based MDR (OR- 
MDRv.1.2). 23 

Interaction information, and interaction graph 

Interaction information analysis used for validation of MDR 
results. 20,24 26 Interaction information was measured among 2 
given loci and subject's status (asthma or control) using Shan- 
non entropy. 27 Let H(X) be the Shannon entropy of X. The infor- 
mation gain (IG) was derived as follows: 

IG(ABC)=I(A; B|C) - I(A; B) 

I(A; B|C)=H(A|C) + H(B|C) - H(A; B|C) 

I(A; B) =H(A) + H(B) - H(A; B) 

, where I (A;B) denotes the dependency of correlation between 
A and B, and I(A; B|C) denotes the interaction of A and B given 
C. When the difference between these 2 parameters, IG (ABC), 
is positive, it is defined as synergy or evidence of an attribute in- 
teraction, whereas when interaction information is negative, it 

Table 2. Linkage disequilibrium among 7 SNPs of 3 MCP genes 



is defined as redundant or evidence of an independency. 20,24,25 
An interaction graph was drawn up based on interactions infor- 
mation with software Orange. 28 

RESULTS 

Hardy-Weinberg cquilbrium, linkage analyses 

The genotype distributions of all 7 single nucleotide polymor- 
phisms (SNPs) were in Hardy-Weinberg equilibrium (data not 
shown). MCP1-25180A had complete linkage (|D'| = 1.0) with 
MCP1 -2076 A > T and with three SNPs {MCP3-420T> C, MCP3+ 
136C> G, and MCP3+5630 T) of the MCP-3 gene in both asth- 
ma and normal control. The remaining SNP, MCP3+984G>A, 
also showed strong linkage with MCPi-25i8G> A MCP1-2076A> 
7/also had complete linkage with MCP2+46A > C, MCP3+5630 T, 
and MCP3+9840A in both asthma and normal control (Table 2). 

Haplotypes distribution and association analyses 

We constructed haplotypes with 2 SNPs of MCP-1 gene, and 4 
SNPs of MCP-3 genes. Only 3 haplotypes were existed among 4 
theoretical haplotypes in MCP-1 gene. Haplotype GA and AA 
comprised 94.0% of all haplotype in MCP-1 gene. Only 5 haplo- 
types were also found among 16 theoretical possible haplotypes 
in MCP-3 gene (Table 3). Haplotype TCCG comprised 82.7% of 
all haplotypes in MCP-3 gene and the each proportion of the 
other haplotypes was less than 10% (Table 3). There was no 
SNP that showed significant association with asthma. There 
was no haplotype that associated with asthma, either (Table 4). 
We tried to make haplotypes simultaneously with all 7 SNPs in 
3 MCP genes because 3 MCP genes are located in same chro- 
mosome 17, then 12 haplotypes were constructed (Table 3). 
The haplotypes that contains TCCG of MCP-3 gene haplotype 
(GA-A-TCCG, AA-C-TCCG, GA-C-TCCG, and AA-A-TCCG) 
comprised 81.1% of all haplotypes. There was no haplotype 
among 12 that significantly associated with asthma (Table 5). 

Gene - gene interaction models and information analyses 

The combination that composed of MCP1-25180A, MCP2+ 
46A>C, and MCP3+5630 T was the best model to predict the 
asthma (accuracy =0.58, CVC= 10/10, P=0.035, OR=1.96; Ta- 





MCP1-2518G>A 


MCP1-2076A> T 


MCP2+46A>C 


MCP3-420T>C 


MCP3+WOG 


MCP3+5630 T 


MCP3+984G>A 


MCP1-2518G>A 




-j -X- 


0.12* 


1* 


1* 


1* 


0.89* 


MCP1-207BA> T 


1* 




1* 


0.07 


0.04 


1* 


1* 


MCP2+46A>C 


0.12* 


1* 




1* 


0.79* 


1* 


1* 


MCP3-420T>C 


1* 


1 


1* 




1 


1 


0.09 


MCP3+1360G 


1* 


0.12* 


1* 


1 




1 


0.03 


MCP3+5630 T 


1* 


1* 


1* 


1 


1 




1 


MCP3+984G>A 


1* 


1* 


0.74* 


1 


1 


0.62 





*P<0.05; P value is D'from 0 by Chi-square test; The values of |D'| of asthmatics are in up-right and of normal control are in down-left area. 
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Table 3. The list of haplotypes and frequencies 



G GflG 


Ma nlnt\/no 
nd|JlULy|Jt; 


% 

/o 


MCP-1 


GA 


62.3 




AA 


31.7 




AT 


6 


MCP-3 


TCCG 


82.7 




TCTG 


7.6 




TCCA 


3.9 




TGCG 


3.8 




CCCG 


2.1 


MCP1+MCP2+MPC3 


GA-A-TCCG 


41.56 




AA-C-TCCG 


16.44 




GA-C-TCCG 


15.38 




AA-A-TCCG* 


7.69 




AA-A-TCTG 


7.56 




AT A TPPA* 
A I -A- 1 UbA 


q pi 
J.O 1 




GA-A-TGCG 


3.25 




GA-A-CCCG 


2.00 




AT-A-TCCG* 


1.63 




AT-A-TGCG* 


0.56 




GA-C-TCCA 


0.06 




GA-C-CCCG 


0.06 



Table 4. Association analyses of MCP-1, -2, and -3 genes polymorphisms be- 
tween asthma and normal control 



Haplotypes were inferred using EM-algorithm. 

*Haplotype that bears the allele A of MCP1-2518G>A, A of MCP2+46A>C, 
and C of MCP3+563C> T, those were detected as the best 3-loci model to pre- 
dict asthma. 

ble 6; Fig. 1). The next best model was a 2-loci combination 
composed of MCP1-2518G>A and MCP2+46A>C (accura- 
cy=0.56, CVC=9/10, P=0.043, OR=1.56; Table 6; Fig. 2). Only 
14 genotype combinations were found among 27 theoretically 
possible combinations in the 3-loci best model. High-risk and 
low-risk genotype combinations for asthma and odds ratios 
(OR) of each genotype combination are designated in Fig. 1. 
More high-risk genotype combinations for asthma were ob- 
served in subjects with the CC alleles of MCP3+5630 T, AA or 
AC alleles of MCP2+46A>C, and AA or GG alleles of MCP1- 
2518G>A. Of the 14 combinations, 6 genotype combinations 
had an OR >1.0. The combination of AA-AA-CC (in order of 
MCP1 -251 8G >A, MCP2+46A > C, and MCP3+5630 T) was the 
best genotype combination, and the OR was 3.67 (95% CI, 1.18- 
13.04). The combination of GA-AC-CT (in order of MCP1- 
2518G>A, MCP2+46A>C, and MCP3+5630T) was the sec- 
ond highest risk combination and the OR was 3.0 (95% CI, 0.97- 
11.0; Fig. 1). 

To confirm the influence of strong linkages among some SNPs 
on the process of detecting for gene - gene interactions, addi- 
tional analyses were carried out with the extra-combination of 
haplotypes that bears the allele A of MCP1-2518G>A, A of 
MCP2+46A>C, and C of MCP3+563C> T. This combination 
was composed of 4 haplotypes that constructed with all 7 SNPs 



Gene 


Polymorphism 


Genotype 


Asthma 


Control 


Pvalue 


OR 


MCP-1 


MCP1-2518G>A 


GG 


166 (41.4) 


143 


(35.8) 


0.39 


ref 






GA 


180 (44.9) 


200 


(50.4) 




0.78 






AA 


54(13.7) 


57 


(13.9) 




0.97 




MCP1-2076A> T 


AA 


360 (90.1) 


347 


(86.9) 


0.32 


ref 






AT 


39 (9.7) 


53 


(13.1) 




0.74 






TT 


1 (0.2) 


0 


(0) 




0 




Haplotvpe 

PA 

bA 


/ 

7- 


CC HO 71 
DO ( 1 6.1 J 


53 


(13.3) 


U.ZZ 


TGT 






uA/- 


1 7Q IA A J\ 


201 


(50.4) 




1 1 K 
LID 






PA /OA 

bA/bA 


lb/ [41./] 


145 


(36.3) 




n on 




A A 

AA 


7- 


1 0"7 / A O O \ 

la/ (4a. o) 


187 


(46.7) 


u.ya 


TGT 






AA/ 
AA/- 


1 CO lOQ C\ 

i Do loa.oj 


172 


(43) 




1 1 A 
1 . 1 4 






A A /A A 

AA/AA 


A A /1 1 1 \ 

44 [I 1. 1) 


41 


(10.4) 




n nn 




AT 

Al 


7* 


obU (auj 


347 


(86.7) 


U.ll 


TGT 






AT/ 

Al/™ 


oa (a. I) 


53 


(13.3) 




n 1 a 
U. 14 






AT /AT 
AI/AI 


1 (u.z) 


0 


(0) 




n 
U 


Ivlur-Z 


IVILr/+4o A > L 


A A 

AA 


909 /cn c\ 
ZUZ (DU.bJ 


170 


(42.1) 


n n7 
U.U7 


i6T 






AP 

AL 


icq / /19\ 

IbO 


190 


(47.9) 




n pi 
U.o 1 






PP 
LL 


9D 11 C\ 

oil l/.OJ 


40 


(10) 




n C9 
U.bz 


H/jpp q 


IVILrj-4/U 1 > L 


TT 
1 1 


oo4 (ab) 


387 


(96.7) 


U.Z4 


TGT 






PT 
L 1 


lb (4) 


10 


(2.5) 




I .b I 






PP 

LL 


U (U) 


3 


(0.8) 




n 
U 






PP 


J4 ItS.Dj 


27 


(6.7) 




1 99 






GG 


0(0) 


0 


(0) 




0 




MCP3+5630 T 


CC 


345 (86.4) 


337 


(84.2) 


0.82 


ref 






CT 


52(12.8) 


60 


(15) 




0.84 






TT 


3(0.8) 


3 


(0.8) 




0.97 




MCP3+984G>A 


GG 


374(93.6) 


370 


(92.5) 


0.74 


ref 






AG 


25(6.1) 


30 


(7.5) 




0.81 






AA 


1 (0.3) 


0 


(0) 




n 
U 




Haplotvpe 
TCCG 


-/- 


14(3.7) 


12 


(3.2) 


0.76 


ref 






TCCG/- 


111 (27.7) 


110 


(27.4) 




1.01 






TCCG/TCCG 275(68.6) 


278 


(68.5) 




1.12 




TCTG 


■I- 


373 (93.3) 


367 


(91.8) 


0.61 


ref 






TCTG/- 


25(6.3) 


33 


(8.3) 




1.01 






TCTG/TCTG 


2(0.5) 


0 


(0) 




0 



ref, referent; OR, Odds ratio; Numbers in parenthesis indicate the percentage; 
Symbol of - indicates the absence of the applicable haplotype; Logistic regres- 
sion with adjustment of age and sex was used to calculate Pvalues and odd 
ratios. 

of 3 genes, including AA-A-TCCG, AT-A-TCCG, AT-A-TCCA, and 
AT-A-TGCG (Table 3 and 5). There was no association between 
the extra-combination of haplotypes and asthma (Table 5). 

The interaction information analysis revealed that MCP1- 
2518G>A and MCP2+46A> Chad positive synergy and provid- 
ed more information (0.82% extra information for the subject's 
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status) than its single effect (main effect, 0.52% vs 0.64%, re- 
spectively, MCP1-2518G>A versus MCP2+46A>C) for a sub- 
ject's status (asthma or normal). MCP3+1360G and MCP3+ 
984G>A also showed positive synergy and provided 0.43% ex- 
tra information for the subject's status than its single main ef- 
fect (Fig. 3). More high-risk genotype combinations for asthma 
were observed in cases bearing the AA or GG alleles of MCP1- 
2518G>A or AA or AC alleles of MCP2+46A> C, similar to the 
3-loci model (Fig. 2). The other SNP combinations had lesser 
synergistic effects than their single main effects. 

DISCUSSION 

Many studies suggest that MCPs are involved in asthma, 5 " 11,29 

Table 5. Association analyses of haplotypes that constructed simultaneously 
with all 7 SNPs between asthma and normal control 



Haplotype 




Asthma 


Control 


Pvalue 


Odds ratio 


GA-A-TCCG 


-/- 


119(29.8) 


131 (32.8) 


0.97 


ref 




+/- 


213(53.3) 


222(55.5) 




0.95 




+/+ 


68(16.9) 


47(11.7) 




0.63 


AA-C-TCCG 


-/- 


289(72.2) 


269(67.2) 


0.69 


ref 




+/- 


104(26) 


117(29.3) 




1.21 




+/+ 


7(1.8) 


14(3.5) 




2.12 


GA-C-TCCG 


-/- 


298(74.4) 


282(70.4) 


0.43 


ref 




+/- 


93(23.3) 


101 (25.3) 




1.15 




+/+ 


9(2.3) 


17(4.3) 




1.98 


AA-A-TCCG 


-/- 


331 (82.7) 


354(88.4) 


0.07 


ref 




+/- 


64(16) 


43(10.8) 




0.63 




+/+ 


5(1.3) 


3(0.8) 




0.57 


AA-A-TCTG 


-/- 


346(86.4) 


339(84.7) 


0.89 


ref 




+/- 


51 (12.8) 


58(14.5) 




0.47 




+/+ 


3(0.8) 


3(0.8) 




0.98 


EC 


-/- 


295(73.8) 


300(74.9) 


0.88 


ref 




+/- 


94(23.4) 


97(24.3) 




1.02 




+/ + 


11 (2.8) 


3(0.8) 




0.28 



N (%); ref: referent; Symbol of - means the absence of haplotype and + means 
the presence; EC: Extra-combination includes AA-A-TCCG, AT-A-TCCG, AT-A- 
TCCA, and AT-A-TGCG that bears the allele A of MCP1-2518G>A, A of MCP2+ 
4BA> C, and C of MCP3+563C> T, those were detected as the best 3-loci mod- 
el to predict asthma. 



thus, we attempted to find an association of SNPs of three MCP 
genes with asthma. No single SNP of MCP-1 or MCP-2 gene was 
associated with asthma. Our results accorded those of a previ- 
ous study indicating that no association existed between MCP- 
3 gene SNPs and asthma. 30 Haplotype constructions were per- 
formed with all 7 SNPs of 3 MCP genes, as well as separately, 
because 3 MCP genes are located in same chromosome 17. 
There was no haplotype that associated with asthma in MCP-1 
and MCPS genes either. 

A gene - gene interaction is often described as "epistasis," 
which is defined as a nonindependence of effect. In other 
words, when the effect of one locus is not predictable unless the 
value of another locus is known, these effects are referred to as 
"epistatic." 31 Many of the important biological consequences of 
gene interactions depend on specific locus-to-locus interac- 
tions, rather than on some average level of interaction. 31 There- 
fore, we hypothesized that a single locus can affect another lo- 
cus, even this locus does not show the efficiency on the expres- 
sion of phenotype alone, and attempted to find the gene - gene 
interaction among SNPs of three MCP genes in asthma. 

Gene - gene interactions were tested using 2 different meth- 
ods in this study. We constructed a gene - gene interaction mod- 
el using a MDR method. The combination of three MCP gene 
SNPs {MCP1-25180A, MCP2+46A>C, and MCP3+5630T) 
was found to comprise the best model. Using this model, we 
could determine which genotype combinations contributed 
more to asthma. Of the 14 genotype combinations indicated in 
the best three-loci model, nine genotype combinations were 
classified as high-risk for asthma, and these combinations in- 
cluded the CC genotype of MCP3+5630 T and the AA or AC 
genotype of MCP2+46A > C with MCP1 -251 8G>A more so than 
others. Among nine genotype combinations, the combination 
of AA-AA-CC (in order of MCP1-2518G> A, MCP2+46A> C, and 
MCP3+5630T, OR 3.67) and GA-AC-CT (in order of MCP 1- 
2518G>A, MCP2+46A>C, and MCP3+5630T, OR 3.0) had 
higher odds ratios to asthma (Fig. 1), which suggests that sub- 
jects with these genotype combinations may be at higher risk 
for asthma. 

We analyzed for linkage disequilibrium separately with asth- 
ma and normal control. But there is no difference in linkage 
disequilibrium among 3 SNPs in which was detected as a best 



Table 6. Gene - gene interaction models 



No. of loci 


Combination of SNPs 


Accuracy 


CVC 


Pvalue 


OR (95% CI) 


1 


MCP2+46A>C 


0.54 


6/10 


0.081 


1.41(1.07-1.87) 


2 


MCP1-2518G>A, MCP2+46A> C 


0.56 


9/10 


0.043 


1.56(1.18-2.06) 


3 


MCP1-2518G>A, MCP2+46A> C, MCP3+5830 T 


0.58 


10/10 


0.035 


1.96(1.46-2.62) 


4 


MCP1-2518G>A, MCP1-207BA> I MCP2+46A> C, MCP3+5B30 T 


0.59 


5/10 


0.812 


2.10(1.58-2.81) 


5 


MCP1-2518G>A, MCP1-2076A> I MCP2+4BA> C, MCP3-420T> C, MCP3+5630 T 


0.6 


9/10 


0.071 


2.35(1.77-3.14) 



P values were validated based on 1,000 permutation test. CVC, cross-validation consistency; OR, Odds ratio; CI, Confidence interval. 
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CC 
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OR 1.0 

8 8 


OR 0.69 
22 JL 




OR 0.8 
12 15 


OR 3.0 
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MCP1-2518G>A 
GA 



GG 



OR 1.0 

3 3 



















Accurary0.58, CVC 10/10, P=0.035, OR 1.96 

Fig. 1. The best model composed of MCP1-2518G>A, MCP2+46A> C and MCP3+563C> T; In each cell, the left bar (orange) represents an asthma, and right bar 
(blue) a normal control. High-risk genotype combinations are shaded dark grey, while low-risk are shaded light grey; OR of each genotype combination was designat- 
ed in each cell; P value was validated based on 1 ,000 permutation test. CVC, cross-validation consistency; OR, odds ratio. 



MCP1-2518G>A 

AA GA GG 

OR 1.57 I OR 0.99 I OR 1.45 
90 91 87 

- - I I I ■ 



OR 1.13 OR 0.73 105 OR 1.1 




Accurary0.56, CVC 9/10, P=0.043, OR 1.56 



Fig. 2. The 2-loci model composed of MCP1-2518G>A and MCP2+4BA>C; In 
each cell, the left bar (orange) represents an asthma, and right bar (blue) a nor- 
mal control; High-risk genotype combinations are shaded dark grey, while low- 
risk are shaded light grey; OR of each genotype combination was designated in 
each cell; Pvalue was validated based on 1,000 permutation test. CVC, cross- 
validation consistency; OR, odds ratio. 

model between asthma and control. These 3 SNPs showed 
strong linkage (|D'| = l.O) in both asthma and control. This sug- 
gests that linkage disequilibria did not affect the MDR analysis. 
In addition to this, to confirm whether these strong linkage dis- 
equilibria were influence on the process of detection of epistat- 
ic models, we made the extra-combination with 4 haplotypes 
(AA-A-TCCG, AT-A-TCCG, AT-A-TCCA, and AT-A-TGCG) that 
bears the allele A of MCP1-2518G>A, A of MCP2+46A>C, and 
C of MCP3+5630T. No significant association between the ex- 




Fig. 3. Interaction graph among seven SNPs of MCP-1, -2, and -3 genes; Boxes 
indicate the main effect of each locus on subject's status, and connections indi- 
cate the interactions effect between two loci on subject's status; Main effect 
means individual attribute to asthma-control status; Interaction effect means 
pairwise combination of attributes to asthma-control status; Solid lines indicate 
synergistic effect between two SNPs to asthma-control status. Dotted lines in- 
dicate independence effect of two SNPs to asthma-control status. 
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tra-combination and asthma indicates that linkage disequilib- 
ria did not affect the process of detecting gene - gene interac- 
tions. 

To validate the MDR results, we used the measure of interac- 
tion information to provide a statistical interpretation of the 
gene - gene interaction models. 20,24,25 Two combinations 
{MCP1 -2518G>A, MCP2+46A > C and MCP3+1360 G, MCP3+ 
984G>A) showed a greater positive synergistic effect on a sub- 
ject's status (asthma or normal) than an individual contribu- 
tion. Interestingly, the combination of MCP1-2518G>A and 
MCP2+46A>C was already determined to be the best two-loci 
model in the MDR analysis, although it was not the best model 
(Table 4). The high-risk genotype combination of this 2-loci 
model was similar to that in the best three-loci model. The 
combination of AA in MCP1 -251 8G>A and AA or AC in MCP2+ 
46A> C, or the combination of GG in MCP1-2518G>A and AA 
or AC in MCP2+46A> C composed the high-risk asthma group, 
with an OR > 1.0 (Fig. 2). MCP1-2518G>A was reported associ- 
ated with asthma. 11,32 Furthermore, it was also reported to be 
associated with various lands of inflammatory disease such as 
bronchiolitis obliterans, sarcoidosis, HIV infection, and metabol- 
ic syndrome in many races. 33 36 However, the MCP1-2518G>A 
polymorphism is not a risk factor for near-fatal asthma and 
shows genetic heterogeneity in different races with respect to 
regulation oiMCP-1 expression. 32 Nevertheless, MCP1-2518G>A 
is still one of the important candidate polymorphisms, and 
MCP genes are important in the development of asthma. 

We cannot conclude which model is better between two-loci 
{MCP1-2518G>A and MCP2+46A>C) and three-loci {MCP1- 
2518G>A, MCP2+46A>C, and MCP3+5630T). Regardless 
which model is best, in terms of epistasis, 31 the effect of some 
SNPs of MCP genes may not be disclosed if the effects of other 
SNPs are not considered. MCP1-2518G>A and MCP2+46A>C 
combination gave us most the information regarding a patient's 
status, and this combination was included in the best three-loci 
model. We could not measure the attributable information gain 
(IG) among three loci with the current method. Moreover, nei- 
ther an appropriate method to reveal gene - gene interactions 
among three different loci nor a knockout animal model of sev- 
eral MCP gene loci exists to validate our findings biologically. 

In summary, gene - gene interactions among SNPs of 3 MCP 
genes in asthma revealed novel findings. Our findings on the 
epistatic influence of the SNPs of three MCP genes suggest that 
these variants may be possible markers for asthma. Further- 
more, they provide important evidence for understanding the 
genetic background of allergic diseases and may aid in the de- 
velopment of a possible genetic marker to predict the develop- 
ment of asthma. 
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